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ABSTRACT 

In this thesis entitled ’COMPUTER AIDED PARAMEtER 
EXTRACTION OF MOSFETS* , the parameters characterising the 
electrical behaviour of N and P channel enhancement mode 
MOSFETS have been determined* 

The device model chosen is the one developed by J*E. 
Meyer [l]. The model when extended to the device under study 
has three parameters and a single independent variable. 

The data points were obtained from the experimentally 
observed I-V characteristics of the transistors. 

The parameters were extracted using the Nonlinear least 
square curve fitting algorithm suggested by Levenberg Marquardt 
[2], The program for the same was developed as a part of the 
project. 

The parameters thus extracted were put to rigorous 
validity checks including checks with l/f noise measurements 
to ascertain their physical consistency. The results justify 
the scheme adopted for study* 
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CHAPTER 1 


INTRCDUCTION 

The CMOS technology i* here to stay, as it finds 
tremendous favour with the manufacturers of LSI and VLSI 
components. It is not surprising therefore that considerable 
work is being done in this field to study the device 
behaviour over a myriad set of operating ranges and variables 
Many aspect of device behaviour continue to elude theore- 
tical explanations. Various models of device have been 
suggested and vdll continue to be formulated to define the 
device behaviour as accurately as possible. 

The behaviour of the device is determined by the 
external variables such as temperature, voltage, current, 
circuit configuration in which the device is used, type of 
load etc., as also by the physical parameters of the device, 
such as doping density, geometry of the device, type of 

semicenductoxs used etc. while the external variables are 

Is 

easily determinedjit 4'the physical parameters of the device 
which are not so easily found. The manufacturers either 
by design ox accident do not dwell on the physical parameters 
that would help dsfine the device behaviour in totality. 
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The converse is also true. If a device with a 
certain behaviour be required, it is absolutely essential 
to establish a cause-effect correlation during the 
manufacturing process . The ’cause’ corresponding to the 
device parameters and ’effect’ to the device behaviour. 

Therefore schemes have to be developed by users and 
manufacturers of the device to extract the device parameters, 
and thuj provide a bridge over a vital information gap. 

Breadboarding used to be a common scheme for parameter 
extraction of discrete devices, since these entailed fewer 
variables, simpler models and usually operated in domains 
whore theories are well established. However the advent of 
L3l and VLSI technology has as a matter of need ushered in a 
new era of computer aided simulation. The models required 
have to he very accurate, since a slight inaccuracy in a 
single stage of the chip, may have a drastic effect on the 
overall behaviour of the chip. The accuracy desired in 
parameter extraction for CAD application too, is of a high 
order. 

The device chosen for the study in this project is 
CMOS inverter GD4007. It has three inverters of which one 
is internally connected and in the remaining two ITMOS and PMOS 
pair are left without internal connection. Some data has also 
been taken on the high speed CMOS 74HC04, 
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The device model chosen in the case of CD4007 is 
the one developed by J.ll. Meyer [ij. The model is discussed 
in Chapter 2, boras basic device physics concepts shaping 
the evolution of this model axe described in Appendix. A. 

The I-V relationship as propounded by J.B. Meyer 
fox N and P. channel are nonlinear within the parameters. 

Hence the parameteri extraction algorithm has to be a 
nonlinear one. The algorithm chosen for parameters extraction 
is the one suggested by Levenberg-lvUrquardt [2j,[3], The 
Levenberg-Maxquardt algorithm is a combination of Gauts-NawWn 
and steepest descent methods. This algorithm provides the 
best tradeoff betv/een computation time and accuracy of 
initial guess to be provided to it as input. The optimisation 
scheme is discussed in Chapter 3. The algorithm for the 
program is listed at Appendix B, 

The parameters thus extracted are discussed in 
Chapter 4, . A comparative study of the expected values and 
obtained values are carried out. The consistency or lack 
of it for the obtained parameters, over a different set of 
data points are also described in this chapter. The validation 
and error analysis are also carried out in this chapter. 

The validation of the parameters have been carried out by 
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a) nuraericelly generating the I-V characteristics of 
N channel and P channel, and ensuring point to 
point coincidence with experimental I-V chara- 
cteristic# 

b) numerically generating the transfer characteristics 
of CMOS inverter comprising of the two transistors 

and comparing it with actual transfer characteristics. 

Also the spectral power density for l/f noise [4] of N 
channel Mosfet has been analysed. This analysis shows a fair 
agreement , with established theories. In Chapter 5 the 
entire work has been summarised and suggestion for further 
work in this field have been incorporated. 



CHAPTER 2 


MODEL OF IviOSFETS 

For computer aided simulation one of the most 
important ingredients is the modal of the device. While 
there are no dirth of models that have been suggested for 
MOSFETS, the choice of a particular model depends mainly on 
the 'aspects of study' of the device one is interested in, 
be it D.C., A.C* transient, noise behaviour or any other 
field. Some models may involve parameters that may cover the 
wide range of operation of the device, while othe:^ have been 
formulated with a certain bias to the field of study. 

Further one will do well, to chose those models which have 
a well established theory rather than a pure empirical model. 

In this project the model chosen for parameter extra- 
ction, is the one developed by J.E. Meyor. In this chapter 
we will focus our attention to some of the salient features 
of this model. The basic device physics concepts contributing 
to the evolution of this model have been listed in Appendix A, 

2.1 TRIODE REGIOM: 

The structure of a typical N channel Mosfet has been 
illustrated in Fig, 1, Let the potentials at gate drain, 
source and substrate be and V|^ respectively, Furt,h€r 



liviiii 
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VD, 
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let V(r) be the p>otential at the surface of the channel and 
Wjl*, and T be the channel width, source drain spacing, 
and ©xide thickness respectiv®ly» ® distance ’Z’ from the 
source the total charge per unit area Q , induced in the 
silicon substrate will consist of a mobile negative charge 
in the surface inversion layer and a fixed negative 
charge Q^, in the underlying depletion region, due to 
ionised impurity atoms* 

Thus 

QjuCx) * Qg(z) - Q^(z) ... ( 1 ) 

By Gauss* iaw 

E 

Q.(z) * - (V -Vtx-.V(z)) (2) 

* ^ox ^ 

and 

Q^{z) * - f 2 ESiqN(V(z)-Vj^+20F) (3) 

Equation (2)Teiat*s the total charge in the silicon to the 
impressed field V -V(z) and an effective threshold voltage 
^tx * T^zs effective threshold voltage includes the 
contribution of metai-semiconductor work-function difference 
Fermi potential of the substrate 0F and the 
equivalent insula ter charge Q^-. Thus, 

''tx = 2 

OX 


(4) 
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£qn, (3) relates the charge of ionised impurity to the 
channel-substrate potential. The voltage drop across an 
elemental section of the channel is given by 


dv ts idR 


V^IaQih(z) 


(5) 


substituting equations (i),(2),(3) and (4) in equation 
(5) and integrating 


L 

- / Id dz 

o 


,/ '* [-(V,-V^,-V(z))) 


ox 


^ 2 -::SiqiM(V(z)-V.+20f]dV 


(6) 


ox 


Id . KC-(Vg-V^^-V(z))* - |R(VCz)-Vj,+20f)®/2j'^d (7) 


where 


K - ItSlS 

■ 2LT„ 


(3) 


f 2E SiqN 


(9) 


-ox 


2.2 SATURATION REGION: 

Eqn. (7) establishes relationship of drain current 
with V »Vi-»V and 1*» triode region of operation. The 

^ O S Q 
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saturation region behaviour demands that the entire channel 
be inverted. This implies that the entire terra in square 
bracket in eqn. (6) (the mobile channel charge) must be 
negative for all z, thus giving the following relationship 

S (io3 


The saturation voltage can be determined by solving 

the above equation by substituting for V and applying 

equality. 


Vsat determined, defines the drain voltage 

at which Qjjj(i) the drain becomes zero. Since in 
practice cannot become zero, implying zero drain 

current, so an approximation has been made. This approxi- 
mation indicates that beyond pinchoff, the drain conductance 
becomes zero while the drain current remains constant at 
a value equal to that achieved at pinch off. However as 
will be noted from typical device characteristics, the drain 
conductance is by no means zero in the saturation region. 
This conductance is generally attributed to two feedback 
mechanisms by which a change in drain voltage can effect the 
drain current* 

The first mechanism is due to a spreading of 
depletion, region near drain, which results in a reduction 
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of active channel length* The second mechanism by which 
a change in drain voltage can effect the drain current 
is by electrostatic coupling of the drain contact region 
to channel. For the purpose of computer aided analysis it 
is usually sufficient to define the drain current empirically. 

Ids = ( 1 + Y . V35) (U) 

When Iq is drain current at Vq * ^sat* 

2.3 MOBILITY REDUCTION: 

With the increase in gate source bias, or . drain 

is a decrease in inversion layer 

mobility due to diffuse surface scattering. This reduction 

in 

in mobility can be accounted for/an empirical fashion. 

^eff =* 

s 

4 sT 

when E = 6x10 V/cm 

SO 

p^^^(o) SI Constant (low field) mobility for 

q « an empirical constant (usually between ,08- ,25) 

Vnc 

E * -(V -V. - )/T (13) 

s ' QS tx 2 ox ^ ^ 

2.4 CURRENT A^D VOLTAGE RELATIONSHIP OF MEYER'S MODEL WHEN 

EXTEiOED TO DEVICE UlOER STUDY: 

The J.E, Meyer's model also account for voltage 
dependent capacitances. Since these are not pertinent to the 
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extraction of parameters determining the d,c, behaviour of the 
device, no attempt is made here to dwell on these. 

The device used for study is CD-4007. On this I,C. 
all substrates and source for both N and P channel are 
internally connected. To determine the I-V characteristics 
of the device these in turn have been grounded. Thus we 
have 

V, . Vj, = 0 (14) 

Solving eqn. (7) we have for the N channel 

- I R(Vi3s+20f)^/^3 Ci5a) 

The drain current for the P channel 

ho - K =[ fR(2(!Sf)3/2 . Vso^t- 2 (V 52 -V^^)V^O 

- 3 (Vj_[,+2«lf)^'^^] (I5b) 

The saturation voltage for the two N and P channel transistors 
are obtained by solving eqn* (10). 

VDSAT^f^Chl ) = |- - I lf(R2.4V^^-HtVQ3+8(!lf) 

(16a) 



12 


''D&AT<‘’ehJ> • C''SG-Vtx>+ ^ - I ir R^-'Vt,H.4Vs^+3(«f) 

(16b) 

The saturation current for the two are exactly as 
indicated in eqn. (11) i,e. 

= Id (17a) 

IsD^Pchl) = ^daat U+Y.Vso) {17b) 

The mobility degration for both the transistors are 

accounted fox by eqn. (12), which is reproduced below. 

t‘off = W°) 

S 

Thus eqns, (15), (16), (17) and (18) define the d.c. behaviour 
of the two N and P MOSFETS. 


Examining the above equations, we have the following 
parameters, which when determined, can be used to simulate 
the device behaviour. 


K 


u E W 
' QX 


2LT 


ox 


R « iTSqEb-iN 

^OX 


20f 




N. 


In 


^i 


(19) 

( 20 ) 


( 21 ) 



< 


i3 

Y « conductance parameter in saturation 
q = mobility degradation factor 

Vtq * Threshold voltage « R]f20f . (22) 

2.5 EXTENSION OF EXTRACTED PARAIviHTERS TO SPICE: 

In the simulation program SPICE, the MOSFET model 
is derived from the FRQHMAN-GPiOVE model. The d.c, 
characteristics of the MOSFET are defined by the parameters, 
VTO, BETA, LAMBDA, PHI and GAMm, where 

GAiviMA « Bulk threshold parameter 

VTO * Zero bias threshold voltage 

BETA *= Intrinsic trans conductance parameters 

PHI » Surface potential at strong inversion 

LAMBDA = Channel length modulation parameter. 

Interestingly J.E,Me'yej.'*s model is used in SPICE 
to determine the voltage dependent capacitance. 

If we examine the inputs to SPICE we find that these 
are either exactly identical ox at best have a minor 
variation to the parameters being highlighted in Meyer’s 
model, (For the simulation of a device, the above mentioned 
parameters have to be provided to SPICE. This is not easy, 
and hence the SPICE assumes certain default values, which are 
bound to be in error with the actual device parameters and 
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the device eimulation at best is approximate)* A quick 
look at the table below establishes the correspondence 
between the SPICE model parameters for MOi^FETS and the 
MOSFET parameters highlighted by Meyer's model* 


SPICE PAEAJvCTERS RELATION 'WITH MEYER’S 

PARAMETERS 


1) VTO » ZEPQ SIAS TMSHOXD ^ ^ 20f 

VOLTAGE ’ 


2) BETAa INTRINSIC TRANS 

CONDUCTANCE ■ = f(K) where 

PARAlvtETER n p W’ 

^ qx 
K. :» 2LT 

3) PHI» SOUFiCE POTEOTIAL AT = 0 = (PHI) 

STRONG INVERSION 


4. GAMA=* BULK THRESHOLD 
PARAI\4ETER 


f(h) where 


R 



V2q £ s^N 


5) LAMBDA= CaANI'JEL LENGTH 
MODULATION 


accounted for by mobility 
reduction expressions. 



CHAPTEa 3 


nonlinear least SQUAxRH method 

Having provided the model of the device and defined 
the parameters to be extracted* it is imperative to think 
of how to extract these parameters. 'vVhile this problem 
can be attacked from many a direction, the best choice 
probably is to see the problem as one of ’curve fitting’. 

This class of problem offers a multiplicity of choice by 
way of standard algorithm, and a few library subroutines. 

Before one chose* a particular algorithm the compatability 
of the model to the extraction algorithm has to be 
established. 

If we examine equation( 15a) of Chapter 2 we have 

Ij, - K[4/3R(2|if)®/2+2(V2s-Vt^)Vo5-VQg2 

^ 4/3 R(VQ3+20f)^/^] (1) 

This equation is of a classic mathematical form. 

y a F(Xj^,X 2 . A^,A 2 ••• Aj^) (2) 

having K independent variable, and M parameters, A point 

of interest that is revealed by examining equation (1) is 

is 

that the equatien ^nonlinear in parameters (also in independent 



and 

v»xtabie«) >^henc« one will have to choose essentially a 
nofilinear method. The nonlinear methods are iterative in 
nature in which the starting values are picked up, upgraded 
by the algorithm untill a convergence criteria is satisfied. 
Some of the common nonlinear methods used are 

1. Steepest i;;>escent 

2. Geust#Newton algorithm 

3. Levenberg-iviarquandt' s algorithm 

4. Davidson Fletcher and Powell^ s algorithm 
and others. 


Examining these algorithms it was realised that the 
Levenberg-Marquardt* s algorithm offerred the best choice. 

The Levenberg-Marquardt algorithm is a combination of 

, Steepest ^algorithm 

Gaus-Newton and ^ descent, The Leven-berg-Maxquardt’ s £ 

aims to capitalise on the best features of these algorithm 

while rejecting their individual pit falls. 


3,2 THE LEVENBERG-MARQUARDT <S ALGORITHM: 


Regardless of whether the equation is linear or 
nonlinear in the parameters, usually the criteria for 
obtaining the best model parameters are commonly satisfied 
by the least squares objective function. 
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whtn « predicted pr calculated value of dependent 
variable for the ith observation 
N » No. of data points 

» Experimental value of the dependent variable for 
the ith observation. 

The best model parameters axe obtained when the* bbjactive 
function is minimised 

T**'* ** Of J * l»2f..«M (4) 

°AJ 

All algorithms aim to minimise the least square 
objective function and differ only in the route they follow 
to upgrade parameters from one iteration to another iteration. 

Focussing attention to the problems in hand v/e have 
in essence a classical mathematical function, 

Y » F(^,X) (6) 

where » (Xj^,X 2 * ..Xj^) 

Aim is to determine J., that would provide the best fit 
to the experimental data. The Gauss Newton algorithm and the 
Steepest Descent algorithm differ in the rules for 
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upgrading the vector A from kth iteration to (k+l)th 
iteration# The incxeraentai vector 3 A that updates the 
parameter values from one iteration to another is related 
as follows. 

K B ,4A * J^(Y~Y ) (5) 

Zn case of ateapest dMctnt 

B is an identity matrix 
K is a step length* 

This is obtained by assuming only 
the first order terms in the Taylor series expansion 
Y ® f{X,A) at some value of X and A and substituting the 
values in eqn, (4)* 

In case of Gauss Newton 


K » 1 




B « /j 




This is obtained 

by selecting 

the first and second 

order terms 

of Taylor series 

expansion for 

Y - where J 

» Jacobi 

matrix 






• • • •: ■’ 


bA^ 

bk^ 

M 


J » 

< 



(6) 

dA^ 
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J * Transpose of Jacobi 


(Y-Y*) 


- Y 


Y^ - Y, 




N 


(7) 




^ Aj 

Aj - K 

e s 

• s 

-♦ t 


A., - A, 


*M 


‘M 


(3) 


The Levenberg-Marquardt's algorithm combines the two 
algorithms to yield 


(J^J + XD^) A = j"^ ( Y - Y *) 

where..' = Marquardfs parameter 

2 T 

D s scaling vector (J + 0 

^ ~ 1*0 is Nash parameter} introduced to circumvent 

failures when one of the diagonal elements of 
T 

J J is zero. 

^steepest 

In general a ^ , descent procedure is expected 

to converge fox popjc starting values, but requires lengthy 
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»oXution time^ b#c»u5« of itf tendency to cJr4#*»“cro#* 
a valley on the function S, instead of following the floor 
of the valley to a minimurnt The primary oauto of this 
it, that the search directions generated are not linearly 
independent. 

Gaus**»Newton on the other hand will converge 
rapidly for good atarti-AO guesses. In the Lovenberg-Marquardt* s 
algorithm initial value of lambda (\) are large and decrease 
towards zero as the optimum is approached. The rules for 
calculating have been heuristically determined. These are 

a) \ Initial » 001 

b) If the objective function at the end of a given 
iteration; calculated with the upgraded value of 
parameters (S new) is greater than the objective 
function at the beginning of the in iteration 

(ij old) then 

\new » \old ♦ 10.0 

This implies that the upgraded values of parameters at the 
end of a given iteration shows a negative trend to the one 
desiredVsince the least square error is enhanced. By 
increasing lambda (X) one is tending more towards stoopast 


descent. 
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c) If the objective function at the end of the given 
iteration (S new) 1* less than S old then 

Anew *Aold ^ 

This decision reflects the fact that the upgraded para- 
meters show a favourable trend towards convergence. By 
decreasing \ (lambda) one is tending more towards Gauss Newton, 

3.2 SUITABILITY OF THE ALGORITHIvi: 

In 3 typical extraction procedure, parameters are 
usually extracted sequentially, one at a time in small 
groups [5]^ The value of each parameter is assumed fixed 
and accurate for use in extracting further parameters, A 
small portion of the model and data from a limited part of 
devices operating ranges are used in extracting each 
parameters, using the least square minimisation criteria, 

Several difficulties may be encountered with this 
procedure. First the sequential approach to parameter 
extraction does not take into account the interaction of 
parameters. Also the parameters obtained from the limited 
range of the device may not correctly define the device 
behaviour in other range of operation. Another problem is 
that some parameters which have a large impact on the 
simulated characteristics, may be deeply imbedded in the 
model equations. 
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The result of these problem* axe that the model 
does not give the full accuracy of which it is capable of* 

A final difficulty is that the parameters extraction 
routines must be tailored to a specific device model Which 
can produce a significant overhead in modifying or 
replacing a model. 

The Levenberg-Marquardt algorithm is a multi~ 
variable regression algorithm and is thus able to attack 
all the parameters simultaneously. The data points 
selected are for the entire triode region, and thua iho 
usual problems found with other techniques are negotiated 
at the same time. The choice of Levenbexg-Marquardt ' s 
algorithm is thus vindicated, 

3.3 problems peculiar TO MOSFET PAtiAviETEFv EXTIRACTION: 

Most models are based on physical theory but tfeaxe 

are always some parameters which do not have physically well 

others 

defined values, and ^ for which the physical values do not 
yield host fit to the device characteristics, Model parameters 
extraction using experimentally observed data points is 
the obvious solution to tht# probid#* 

In order to be useful oh a routine basis any 
parameter extraction algorithm mutt produce parameters 
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reiiably and accurately. Two potential sources of error 
in the optimisation technique described, are [5], 

a) local minima 

b) Redundancy/insensitivity of parameters 

The first can result in non~optiraum parameter value-if, 
and the second can produce non-unique parameter values, but 
yet give an optimum fit to the data. The latter is 
essentially a feature of the model equation. 

The accuracy in locating the minima is primarily 
determined by the convergence criteria. More strict the 
convergence tolerance, more accurate are the results, but 
this will need more computation time. 

The redundancy problem can be tackled by good initial 
guesses for the parameters; coupled with constrained 
minimisation. Constrained minimisation implies that a b<fund has 
to be placed on upper and lower limits of the parau^ter 
values and also the sign pf the values, 

3.4 THE COiWERGEI'CE CRITERIA s 

In the preceding paragraph we have assextted 
the importance of a judicious choice of convergence criteria.’ 
Some of the common convergence criteria that are usually 
chosen are; 
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a) The objective function value falls below a 
pjedetermined value* 

b) The xelative chajjge in objective function value from 
one iteiation to another fis less than same 
specified value* 

c) The change in each parameters from one iteration 
to the next is within some specified value. 

d) No adjustment of can be found that reduces the 
error. 


The convergence criteria chosen in this case is a combi- 
nation of criteria mentioned in subpaxa (a) and (c) above, i.e. 

- “ < lE-7 

S < la-7 


The above twin conditions have to be jointly satisfied for 
convergence to achieved. The justification for the choice 
is listed below. 


The scheme has been developed for extracting three 


parameters. 


Ml) 


K * 


E 


PI. 


21 T 


■W 


ox 


T 

A(2) R =* YSqESjN 

OJC 

A(3) » 2l5f 



25 


The expected values cf: 

A(l) is in the neighbourhood of iE*.3 
A(2) is in the neighbourhood of J£0 
A(3) is in the neighbourhood of i£0 

Therefore a change in these parameters by a quantity less 
than lH»7 is indicative of the fact that the change is 
less than .001/-, meaning thereby that further refinements, 
except for aesthetics have little influence on the 
parameters. 

Further the number of data points chosen for 

be 

extracting these parameters were taken to/lO. The value 
of dependent variable (^□ 3 /I 5 q) showed a spread from 
(«l'"iO)H-3 amps. An objective function less than IE— 7 
implies an error of less than .1 ma between the data points 
numerically generated and data points experimentally observed. 
This 03^^®^ of error can easily be accounted to the measure- 
ment technique used. 

The two criteria put together are expected to give 
fairly consistent results over various sets of data points. 

Also to check the otherwise divergent tendency of 
the program a constraint [6j was impo-sssd on the sign of the 
parameters. All parameters values extracted were expected 
to be positive. Hence # during the course of a iteration 
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if a i>axaineter value waa found to be neciative it waa reset 
to corresponding positive value before the start of 
subsequent iteration. 

3.5 THE KiOGRAfvl TO BiPLEiViSNT LEVE?®ERG«iViAHQUARDT » S 
ALGORITHM: 

The program realising the chosen algorithm was 
developed as a part of the project. In order to be able to 
use the program on a routine basis, the program was made 
modular in nature. This ensures little overhead in changing 
the models and the program is adequate for parameter 
extraction for the entire class of multivariable nonlinear 
equations. The changes called for are minimal. These axe* 

i) The dimension statement to be appropriately changed 

ii) SUBROUTIi'lE FUNC to be chaBged to incorporate the 
new model, 

iii) bUBFiOUTINH DHRIV ■» to be changed appropriately to 

dY. 

analytical generate / J » 1, ... M 

* i * 1. ... N 

The flow chart of the program is given in Fig. 2, while the 
detailed algorithm is listed at Appendix B. 

The results obtained and some of the peculiarities 
of the program are elaborated in the succeeding chapters. 



CHAPTER 4 


EXPERIJvENTAL REisULT AiO DISCUSSION 

For a parameter extraction scheme like the one being 
followed there are three maj,o.r requirements* These are 

a) Model 

t>) Extraction algorithm 

c) Experimental data 

While requirements listed in (a) and (b) above have 
been di$cu®se<i in the preceeding chapters, we will now 
examine the shape, and form of experimental data to be 
provided as the’ third and final input to thf ^pheme. 

The triode region cuxrent-voltage re|,iptionship 
of N channel MOSFET (identical for P channel! is given by 
Eqn. (15), Chapter 2 which is rewritten here for convenience, 

( 1 ) 

In this equation, as it has been formed, the drain 
current is obviously the dependent variable. The 
independent variables are the gate to source voltage ^GS 
and the drain to source voltage (It may however be 
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noted that the independent variable may change by rearranging 
the torm$ in a different form). 

The LeVcn^®^9“'i'''‘3rquardt algorithm can take care 
of all the three variables and Vg^) at the same time 

However it was decided to obtain data for vrs V^g for 

constant Vqs because of the following reasons. 

a) The Vj^g plot clearly exhibits the two distinct 

phase of operation; the triode region and the 
saturation region. 


b) It is easy to graphically ascertain the pinch off 

voltage once the parameters have been extracted^ 

the accuracy of the parameters can be ascertained by 
calculating as per eqn. (16a) and (b) chapter 2. 

c) It is easy to generate graphically the transfer 
characteristics of a CMOS, by having the I^gvrs v^g 
and IgQ VIS Vgp plot of the constituent N and P channel 

Mosfets. This transfer characteristics then in turn can be 
compared with numerically generated transfer characteristic 
(using extracted values of parameters) for point to point 
coincidence. 


d) By keeping ^GS constant for a set of data the whole 
equation reduces to a simple one dependent and one independent 
variable. The analysis of the equation becomes easier. 
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4.1 EXPERUVEISTTAL AKRAi'^ElvENr Am RESULTS: 

® ) N Channel enhancement Mode Mi^SF^ : 

The scheme employed to determine the vxs ^DS 
data for the N channel MOSFHT is shown in Fig. 3, A fixed 
bias was applied to the gate using a 12 bit DAG, The drain 
supply was connected through a current limiting resistor. 

The source was connected to the virtual ground of the opamp. The 
drain current was measured using an opamp as a CCVS. 

The voltage is a measure of Ipg appropriately scaled. 

The voltage Vp^ and were measured using a 6~ bit 

digital multimeter. The data points thus determined are 
tabulated in Table 2, 

b) P Channel ; 

The scheme for determining the data points of P channel 

is shown in Fig. 4 and is identical to the one used for 

N channel. Here too the source and substrate are internally 

connected and these in turn are connected to ^ the virtual 

negative 

ground of the opamp. A ^ voltage is applied to both 
source and drain. The data obtained is tabulated in 
Table 3. 

The data that has been tabulated have been restricted 
to triode region of operation of both N and P channel 
Mosfets. Further ae may be apparent from the two tables, 



Table 2 


SN 

^GS(v) 

^DS(v) 

^DS(ma) 

SIM 

^GS(v) 

^DS(v) 

^DS(nia) 

1 

4.2045 

.1145 

.1867 

1 

6 .5015 

.114 

.3198 

2 

“*<10"“ 

.224 

.3466 

2 

-do- 

.118 

.3297 

3 

-do*- 

.322 

.4745 

3 

-do- 

.377 

.9932 

4 

-do- 

.439 

.6098 

4 

-do- 

.438 

1.1382 

5 

-do- 

.527 

.6982 

5 

-do- 

.588 

1.4818 

6 

-do- 

.637 

.7925 

6 

-do- 

,736 

1.7967 

7 

-do- 

.777 

.8879 

7 

-do- 

.843 

2.0082 

8 

-do- 

,814 

.9097 

8 

— do — 

1.010 

2,314 

9 

-do~ 

.891 

.9465 

9 

-do- 

1.158 

2.5585 

IQ™ 

_ -do- 

1.020 

.9912 

J:?, 

-do- 

1.326 

2,8087 

1 

5.007 

.0558 

.1211 

1 

7.007 

.126 

.3764 

2 

-do- 

.112 

.2363 

2 

-do- 

.266 

.7697 

3 

-do- 

.241 

.4836 

3 

-do- 

.492 

1.3686 

4 

-do- 

.396 

.7538 

: ^ 

-do- 

.640 

1.7301 

5 

-do- 

.480 

.8869 

! 5 

-do- 

.837 

2.180 

6 

-do- 

.510 

.9316 i 

i 6 

-do- 

1.037 

2.5932 

7 

— do— ■ 

.615 

1.0816 

i 7 

f 

-do~ 

1 .134 

2.8713 

8 

-do- 

.730 

1.2295 

8 

-do- 

1.380 

3.206 

9 

-do- 

*807 

1.3189 

9 

-do- 

1.597 

3.527 

10 

-do- 

.933 

1.4461 

10 

1 

-do- 

1.811 

3.795 

1 

6.007 

.1364 

.3486 

[l 

8.0026 

.127 

.4221 

2 

-do- 

.206 

.5104 

1 2 

-do- 

.235 

,7697 

3 

-do- 

.397 

.9514 

! 3 

-do- 

.519 

1.6258 

4 

-do- 

.416 

.9932 i 

1 ^ 

-do- 

.698 

2.1314 

5 

-do- 

.570 

1.3090 

! 5 

k 

-do- 

.910 

2.6915 

6 

-do- 

.632 

1.4292 i 

! 6 

-do- 

1.106 

3.1723 

7 

— do— 

♦ 766 

1.6735 j 7 

-do- 

1.353 

3.7255 

8 

-do- 

.850 

1.8165 j 8 
1.9993 ! 9 

-do- 

1.513 

4.0483 

9 

-do- 

.965 

-do- 

1.755 

4.4912 

10 

-do- 

1.091 

2.180 

10 

— do- 

2.136 

5.068 
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Table 3 


SN 

^SG( v) 

SD( v) 

^SD(ma) 


^SG( v) 

^SD( v) 

^SD(ma) 

1 

4,003 

.193 

.2556 

1 

6.0067 

.124 

.255 

2 

-do- 

.222 

.2923 

2 

-do- 

.172 

.349 

3 

-do- 

.464 

.5747 

3 

-do- 

.253 

.5095 

4 

-do— 

.645 

.7620 

4 

-do- 

.373 

.744 

5 

-do- 

.756 

.8651 

5 

-do- 

.495 

.974 

6 

-do-^ 

.934 

1.0148 

6 

-do- 

.617 

1.199 

7 

—do— 

1.188 

1.1931 

7 

-do- 

.650 

1.256 

8 

-do- 

1.204 

1.2021 

8 

-do- 

.823 

1.562 

9 

-do- 

1.252 

1.2358 

9 

-do- 

.943 

1.765 

10 

-do- 

1.362 

1.2913 

10 

-do- 

1.030 

1.910 

1 

5.0030 

.1392 

.247 

r 

"l>.50i5” 

.106 

.230 

2 

-do- 

.236 

.407 

2 

-do- 

.206 

.4409 

3 

-do- 

.366 

.628 i 

3 

3 

-do- 

.309 

.6604 

4 

-do- 

,478 

.806 j 

4 

-do- 

.447 

.9365 

5 

-do- 

.557 

.928 j 

5 

-do- 

.623 

1.2862 

6* 

. -do- 

.673 

1.102 ( 

6 

-do- 

.700 

1.429 

7 

-do- 

.756 

1.221 : 

7 

-do- 

.742 

1.5126 

8 

-do- 

.846 

1.347 

8 

-do- 

.888 

1.785 

9 

-do- 

.928 

1.460 

9 

-do- 

.892 

1.814 

10 

-do- 

1.290 

1.913 

10 

-do- 

1.225 

2.381 

1 

5.5031 

.148 

.282 

1 

7.0031 

.223 

.501 

2 

-do- 

.221 

.418 

2 

-do- 

.282 

.626 

3 

-do- 

.398 

.736 

3 

-do- 

.476 

1.044 

4 

-do- 

.462 

.848 

4 

-do- 

.554 

1.205 

5 

-do- 

,570 

1.032 

5 

-do- 

.736 

1.579 

6 

-do- 

.627 

1.128 

6 

-do- 

.913 

1.936 

7 

-do- 

.765 

1.352 

7 

-do- 

1.195 

2.467 

8 

-do- 

,853 

1.4908 

8 

-do- 

1,393 

2.821 

9 

-do- 

.984 

1.690 

9 

-do- 

1.430 

2.888 

10 

-do- 1.444 

2.331 

10 

-do- 

1.462 

2.942 
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the range of been restricted from 4v to 8v to 

obtain the data* The reasons for doing this are listed below. 


i) For lower values of and 4 .u^ ^ ir j 

GS SG, the range of and 

^SD define the triode region of operation for N and 

P channel transistors respectively, are greatly reduced, and 

the data points tend to get very close. Also the order of 

current being measured is less than .2 ma. Hence a small 

relative error in measurement reflects a large percentage 

error in data points. This may lead to erroneous values of 

parameters, when fed to the extraction program. 


ii) For higher values of and the currents 

and IgQ rise above 7.5 ma. This limit exceeds the 
manufacturers specification for both the IC and the opamp 
being used as a CCVS, Further the virtual ground of the 
opamp starts drifting, v^hich will introduce error in 
data points. 


Besides the data points, another input to the 
extraction algorithm are the initial guesses of the 
parameters being extracted. It was noted during program 
execution that for those functions where a unique solution to the 
parameter values were expected, the initial guess could be 
as wild as possible. The solution yielded was well within 
the acceptable error bound (less than .01/^). However as has 
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been pointed out earlier in Chapter 3, the parameter 
extraction for Mosfets suffer from two potential sources 
of error namely, 

a) The local minima 

b) Redundancy/insensitivity of the parameter. 

V^hile the local minima problem was easily resolved, 
the Redundancy problems demanded a fair initial guess. 
Fortunately device specifications were available which helped 
bring the initial guesses close to the expected values. 

4.2 PARAiViETHR EXTRACTION: 

a) For the test model : 

To be able to verify the veracity of the program, 
a mathematical model, somewhat similar to J»£. Meyer ’ s 
model for mosfet was employed. The mathematical model 
had five parameters, and a single independent variable. The 
data points were numerically generated. The output of the 
program is tabulated in Table 4, 

From the table ' it is 

apparent that parameter A(l) is insensitive, and virtually 
remains unchanged. If this parameter’s initial guess was 
changed to ,99-1.01, it showed a tendency to stay there, 
thus yielding fruiteful results. 
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Table 4 


SN 

Para- 

Initial Expected 

Obtained 

Least 

Av. error/ 


meters 

guess 

val. 

val. 

square 

error 

data points 

1. 

A(l) 

a(2) 

•90 

1.0 

.9047 

.472 

2.1E-5 

2. 

2.2 

2,0 

1.981 


3. 

A(3) 

3.4 

3.0 

3.CK336 



4. 

A(4) 

4.8 

4.0 

4.0073 



5. 

A(5) 

5.90 

5.0 

4.9996 




Another observation that was made during the 
parameter extraction of the mathematical model was that if 
the initial guesses were relaxed increasingly, a point was 
reached when the program showed a divergent tendency and 
showed no signs of turning in the correct direction even 
after 200 iterations. 

Also in absence of any constraint on the parameters, 
the parameters a( 1) and a( 2) showed a trend to turn -ve, 
while the error minimisation was still being affected. This 
prompted an imposition of sign constraint in the algorithm, 
whereby if a parameter yielded at the end of an iteration 
was ngative, it was reset to its corresponding positive 
value before entering the next iteration. 
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Tho se lessons were note<i an<J utilised for parauneter 

extraction of Mosfets, 

b) The MQSFET Parameter Extraction ; 


The Moyer’s model focusses attention on three parameters. 
Having appropriately modified the subroutines ’FUNG’ and ’DERIV* 


and substituted the 

test 

data with experimental data, 

» the 

program was executed to yield the following results. 



Table 5 



4,2.b(i); 

oo 

N Channel 




\s 

K 

I NIT EXT 

, R 

INIT 

2 df Least Square 

^ objective 

EXT INIT EXT Function 

Average 

Exrox/Data 

point. 

1. 4.2045 

.309 .337 
E-3 E-3 

1.645 1^645 0.74 0.74 0.6795E-8 

.255E-4 

2. 5.007 

-do- ,332 
E-3 

— do— 

1.645 -do- 0.74 0.234E-8 

.1554E-4 

3. 6.007 

-do- ,310 

E-3 

-do- 

1.645 -do- 0.74 0.2166E-8 

.149E-4 

4. 6.5015 

-do- .299 

E-3 

-do- 

1,645 - - 0.7399 0.396E-8 

.199E-4 

5. 7,007 

—do— .287 
E-3 

— do— 

1.6449 - 0.7398 0,159E-7 

.399E-4 

6, 8.0026 

-do- .270 
E-3 

-do- 

1.6447 - 0.7386 0.2175E-' 

7 .466E-4 
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4*2,b(ii); P Channel Table 6 


VsQ 

1. 4.CX)30 

2. 5.0039 

3. 5.5031 

4. 6.0014 

5. 6.5012 

6. 7,0047 


K R 

I NIT EXT I NIT EXT 


20f 

2ML. 


EXT 


Least sq. Average 
objective Error/Data 
Fiinct ion; . Point . 


189E-3 

.240 0.78 
£-3 

do- 

.240 -do- 
E-3 

do- 

.229 -do- 
, E-3 

do- 

.218 -do- 
E-3 

■do- 

.209 -do- 
£-3 

■do- 

.200 —do— 
E-3 


0.779 0.62 0.6195 
0.7794 -do- 0.6185 
0.7788-do- 0.6170 
0.7795-do- 0,6187 
0.7790-do- 0.6177 
p. 7709-do- 0.5974 


245E-9 

0.497E-5 

190E-8 

0.14E-4 

387E-8 

0.195E-4 

lllE-8 

.llE-4 

226E-8 

.15E-4 

542E-8 

.238E-4 


INIT = EXPECTED VAL/lNITML VALUE 
EXT * EXTPiACTED VALUE 

c) Use of l/f noise for validation: 

The theoretical expression relating the noise spectral 
power density with the device parameter [4] is given by 

, = ^ ^ ( 1 ) 

with K = 
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for weak invenion jpQ^J < and 

^ID _ K_ 

(Cox^+Cit )“ 
and for strong invenion 1 


is - K 

r 2 


N. 


'•D 


f* 


( 3 ) 


4 • O 

indicating that strongly depends on gate bias 

through l/Q„^. 


In the above three expressions 


SfD 5= current noise spectral power density 
Iq =s drain current 

«« 1 -*3 

= trap density (ev .cm ) 

= tunnelling const, for electron 
P ~ q/kT 

k » Boltzmanns constant 
q = electronic charge 
T = absolute temperature 
= channel charge 


To validate the model and associated parameters the following 
steps could be employed. 
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1* Using Meyer’s model (with the values of parameter 
extracted) l/(pQ^)^ can be calculated. This should 

r\ 

have the same slope as , If this is true, value 

KN. lu u 

of — — of eqn. (3) is known. If not the procedure has 
f* 

to be abandoned because may be Vq dependent and 

this has to be experimentally determined first. 

KNf 

2. Use the value of obtained to check the validity 

of eqn. (2), using and as used in Meyer's model 

and ignoring . if the experiment and theory agree, 

the effect of is not important, otherwise methods 

have to be developed to include in the model and 
use some other model than the one suggested by 
J.E, Meyer. 

In the present work step (1) has been checked and the 
results are discussed below. 

Fig, 5 shows the experimental set up used to determine 
Sid/Iq . The results then obtained are tabulated in 
Table 7 and also plotted in Fig. 6. 







V6(v> 
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TablG 7 

Vq = .5V 


SN V. 


Freq- 

uency 


A^/Hz. 


SN 


1 1.867 6 

8 

10 

2 1.8759 6 

8 

10 

3 1.9150 6 

8 

10 

4 1.9911 6 

8 

10 

5 2.1069 6 

8 

10 

6 2.2192 6 

8 

10 


4.845E-14 { 7 

2.688E-14 

1.467E-14 

I 

6.8889E-14 ! 8 
3.448E-14 I 
2.405E-14 j 

1.371E-14 I 9 

8.76E-15 

7.228E-15 

3.96E-15 I 
3.1C9E-15 j 
2.275E-15 j 

i 

I 

3.588E-15 I 
2,359E-15 ; 

1.873E-15 I 

I 

1 

2.204E-18 j 

9.528E-19 

7.83E-19 


V Freq- 
G uency 


2.4169 6 

8 

10 

2.6062 6 

8 

10 

2.9248 6 

8 

10 



1.50E-1S 

7.48E-19 

5.85E-19 

8.157E-19 

4.176E-19 

3.161E-19 

4.87E-19 

1.756E-19 

9.34E-20 
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For strong invention as has been brought out earlier 



A plot of 
Fig. 6. 


^t 

21^ 


rQ 


( 4 ) 


against has been taken and is shown in 


To determine an approximation has been made 

^ = ^1 ( 5 ) 

for small values of V^. 

The term p refers to effective mobility and its value 
has been taken from the value of K extracted, using 
Levenberg-Marquardt » 

^iD 1 

On comparing the curves — ^ vrs Vq and -4^ vrs Vq, 

we find that the two are nearly parallel, and indScating that 

KN. 

is at best weakly dependent on Vq. The constant — j 
thus obtained for various values of are tabulated in 
Table 8. 

In order to check step (2), in weak inversion 
region (order na amps) has to be measured. This could not 
be performed due to the inability of present setup used in 
the experiment. However it was found that the current 
spectral power densities are well defined in subthreshold 
biases too. 
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T 3bl© 8 


f 



(coulomb volt sec) 

6 Hz 

1.875 


.2943x10"^'^ 


1.9150 


,2749x10“^^ 


1.9911 


.6990x10“^^ 


2.1069 


— 3*^ 

.2949x10 


2.2192 


.924xl0~^® 


4.2.2 Some Comments on the Extraction Algorithm Behaviour 
xelevent to Model under Study: 

Some of the common tendencies of the extraction 
algorithm noted are elaborated below. 

a) Susceptability to initial guess: 

Of three parameters K,R and 20f, the parameter 
R and 20f were mostly insensitive during least square 
error minimisation. The values of these parameters obtained 
on convergence showed an affinity to the initial guesses 
provided for these. If the initial guess for these parameters 



46 


was TBlaxed continuously , the -extrsction aloorithm. tended to 
diverge* 

(b) Sensitivity of the parameter : 

Of the three parameter the parameter K was most 
sensitive and showed scant respect for its initial guess. 

The remaining two parameters were in comparision, insensitive 

( c ) Constrained Minimisation ; 

In absence of any constraint on the parameters the 
program often tended to yield negative values for K and R 
while the least square minimisation was still being achieved. 
In order to obtain the correct sign and values of parameters 
a constraint was imposed on the sign of the parameters. This 
helped bring the program execution to order. 

( d ) Choice of Mathematical Function : 

All the above tendencies were found absent, when the 
extraction algorithm was extended to a well defined mathema- 
tical function, where all parameters were independent and 
unique in their values. The extraction algorithm always 
yielded results with an accuracy better than .0!/ and took 
fewer iterations to converge. 

Hence it may be worthwhile to select those models 
where parameters are well defined, independent and unique. 
However this may not always be possible. 
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4.3 THE MOBILITY BEHAVIOUR: 


As exp€cted the parameter K for both N and P channel 
show a decline with increasing value of and as 
expected because the parameter 


K 


I* E 


ox 


2LT 


ox 


W 


The term refers to effective mobility and as anticipated 
it shows a degradation with increasing gate bias. The graph 
showing the mobility behavioux of N and P channel as 
extracted from the parameters yielded by Levenberg-Marquaxdt’ s 
algorithm are plotted in Fig. 7, 

Using equation fox mobility 
•^eff ^eff(o) ^ 


4 / 

where = 6x10 v/cm 
so 

Es = - 0.5 Vds)/T^^ 


an attempt was made to calculate the- value of 'q’ the mobility 

degradation factor. It was observed that value of q thus 

rather 

obtained showed a wide dispersal ^ than being a constant 
as expected. 

'■A'' 
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an effort was made to go in for newer model, 
which could explaxtv -tJitt «ob-©«ivdd mobility b«havi/>ux* It 
was found that the impejioal -rnofiUii ja.»opc-&«<i t>y S.C* Sun and 
James Plummer [7] explained the -observed mobility ^ibcjiomorwMl 
with a fair degree of accuracy* 


The model proposed aims to derive effective mobility 
from d.c, drain conductance gd. The effective mobility 
is related to the drain conductance by (refer Fig. 8 ). 


(L/W)gd 


■^eff ^ q N, 


/ VD 0 


(7) 


inv 


with 


‘I Ninv = °ox V''t - 5 [i-" 


2f2 In N^/N^ 


]*VD 


(B) 


and a = T2 ( ox 


and Lg the Bulk Debye length 



5 n7^ 


(9) 


From the device specifications provided, in conjunction 
with the experimental data, was computed for both N 

and P channel. The results are tabulated in Tables 8 and 9, 


and also plotted in Fig^ 7 
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Table 8 


N Channel 

SN Vq at Vq 

(mho) 

1 4.2045 1. 6305x10" Vy^Ua45 

2 5.007 2.1702x10'"^/Vq=.0558 

3 6.007 2*555x1o“^/Vq=,1364 

4 6.5015 2.8052x10"Vvq=.114 

5 7.007 2. 9873x10“ 126 

6 8.0026 3. 323x10“ Vvp=. 127 


L/W gel 
(mho) 

^Ninv 
at V-.= 
1.85^ 

^eff 

( cxti/V 
sec.) 

326x10”^ 

.0645x10“^ 

505.4 

434x10“'^ 

.0907x10’*^ 

478.5 

SllxlO""^ 

.1165x10“^ 

437.87 

561x10“^ 

.1323x10“^ 

424.0362 

597x10“"^ 

.1467x10“^ 

406.95 

664x10“"^ 

.1760x10“^ 

377.272 


P Channel 


Table 9 


3N 

Vq 

Gd 

(mho) 

Y:hoi 

^e|f 

cmvV sec 

1 

4.003 

1 . 324xl0“^/Vj^= , 193 V 

♦0169x10“^ ,0720x10“^ 

234.72 

2 

5.0030 

1 . 774x 10 “^/Vq= . 1392V 

,02274x10“^ .1028x10"^ 

221.20 

3 

5,5031 

1,905x10"^/Vq=s.148v 

.02442x10“^ .1173x10"^ 

208.184 

4 

6,0067 

2.056xl0“^/Vj^=.124 

.02635x10“^ .1328x10“^ 

198.41 

5 

6.5015 

2.l69xl0“^/Vp=.106 

.02780x10“^ .1478x10'’^ 

188.09 

6 

7.0031 

2 .2466x10“^^^* .223 

.02880x10“^ .1597x10*"^ 

180.338 
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It should be noted that the effective mobility is distijjct 

from the field effect mobility M-pg> which is obtained 

from the transconductance g and is defined by 

m 





( 10 ) 


The relationship between the effective mobility and field 
effect mobility is also shown in Fig* 8. It is apparent that 
the mobility data from g^ and measurements are essentially 
identical at the point of maximum slope (near Vj) . Beyond the 
point of inflection, the field effect mobility is always 
smallert than the effective mobility. Therefore if ppg is used 
in device modelling currents and device switching speeds will 
be under estimated. 

4.4 VALIDATION OF EXTRACTED PARAMETERS: 

The question that readily comes to ones attention is 
how good are the parameters extracted? Are the pararaetera 
extracted capable of simulating the device behaviour in the 
entire region of operation and if so, then with what degree 
of accuracy? Do also the parameters extracted have a 
well defined physical significance. 

To put these doubts to rest, the extracted parameters 
were subjected to a variety of validation schemes. The 
schemes employed are 
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a) simulation, of. dc characteristics for the N and P 
channel MOSEETS and a comparison with -expejimental 
I-V characteristics* 

b) Simulation of transfer characteristics of CMOS 
inverter comprising of the N and P channel MOSFETS 
under study, and comparing it with experimentally 
obtained transfer characteristics. 

4.4(a): Simulation of d,c. Characteristics: 

The I-V characteristics of N and P channel have 
been simulated using extracted value of parameters. The 
two plots, experimental and simulated are illustrated in 
Fig. 9 for N channel and 10 for P channel. They show 
a fair coincidence for both N and P channel* 

It may be argued that such a validation scheme is 
superflours., since the data point fox extraction algorithm 
using least square minimisation have been provided from the 
experimental I-V graph, and that coincidence is only 
expected since convergence criteria has been satisfied. 

As a matter of fact, it is not so. The data points 
for extraction algorithm have been chosen from a limited 
range of device operation namely the triode region. Here 
the simulation is being attempted for both the triode and 
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saturation region. The most critical area of curve fitting 
is the 'Knee’ or pinchoff voltage referred in this 
thesis as 

It was pointed out earlier that in the l-V relation- 
ship proposed by J.E. Meyer, two of the parameters R and 20.f 
were found rather insensitive, despite the fact that least 
square minimisation was being achieved. These two parameters 
tended to track the initial guess. 

However these parameters are dominant in determining 
the ^Knee' , as such there may arise a requirement fox 
upgrading the initial guesses for these parameters to 
obtain the best fit concurrently satisfying the convergence 
criteria of the extraction algorithm. The fact that such 
a curve fitting has been obtained, helps buildup a confidence 
towards the values of the parameters extracted* 

4, 4(b): Simulation of Transfer characteristics of the CNDS 
Comprising of N and P channel Mosfets: ; 

This scheme in essence is a further enhancement of the 
previous scheme. Here the aim is to check the behaviour of 
the two transistor while interacting with each other. The 
flow chart of the algorithm followed is at Appendix C. Also 
refer to Fig, 11 for the configuration in which the two 
MOSFETS are arranged. 
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The simulated txansfer characteristics and experi- 
mental transfer chararcteristics are plotted in Fig. 12. 

Good point to point coincidence is obtained between the two 
plots thereby further validating the parameter values. 

4.4 ERROR ANALYSIS: 

While care has been taken to minimise error at all 
stage, it is not possible to completely eliminate the error. 
There are several potential sources of error and these must 
be identified to decide on the error tolerance of the result. 
These sources of error are 

a) Accuracy of the model 

b) Measurement inaccuracies 

c) Convergence criteria used for extraction algorithm 

d) Instrumentation error 

4.4(a) Accuracy of the Model; 

The triode region behaviour of the NK)SFETS have been 
found to, be adequately well explained by the Meyer’s model. 
However the linearisation of the model in saturation region 
is not accurate and has a inbuilt source of error. The 
model is specially vulnerable at pinch off, when the 
device operation makes a transition from triode to saturation 
region. While the experimental I-V data is smooth over the 
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entire range of operation, the Meyer's model induces an 
abrupt Jump , in current value at pinch off. 

4.4(b): Measurement inaccuracies: 

An opamp based CCVS has been employed. The leakage 
current of the opamp, though of the order of na have been 
ignored. Further the offset voltage tended to increase with 
increasing value of current. This is bound to introduce 
error, which may not be very significant if proper care is 
taken during measurements, 

4.4(c): Convergence criteria used for Extraction Algorithm: 

Fairly stringent convergence criteria has been used in 
this case, but it is by no means zero as is to be expected^ 
This may lead to a very small error. An attempt to further 
minimise the error needs far greater computation time, 

4.4(d): Instrumentation Error: 

The device used for measurement of all voltages is the 
6^ digit multimeter. It was found that if different 
multimeter scales were used to measure the same data, it 
showed a departure which was scale dependent. To minimise 
the error, the entire set of data points were taken on the 


same scale. 



CHAPTER 5 


CONCLUSION 

Computer aided simulation is finding tremendous 
favour with manufacturers and users of LSI and VLSI chips* 

To be able to simulate the device behaviourvaccurate 
device parameters as also the domain of operation (external 
variables) are the major ingredients 

In this project CMOS CD4007 has been studied in d.c, 

operation and related parameters have been extracted. For 

extraction of parameters Levenberg-Marquardt’ s algorithm 

using 

has been developed and inplemented / experimental 1~V 
data points in the triode region. The program developed 
is modular in nature and can be extended to any model 
function having multiple parameter, with minimal changes. 

The values of parameters have been checked by a number 
of schemes. These are: 

a) Simulation of I-V characteristics of N and P 
channel mosfet, 

b) Simulation of CMOS transfer characteristic 
comprising of the two transistors. 


c) Mobility analysis 
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The parameters values obtained show satisfactory 

number 

results as proved by a / of validation schemes employed. 

Error analysis has been made, and also comments on 
the peculiarities of extraction algorithm under certain 
conditions has been described. Some data has also been taken 
on l/f spectral power density of N channel MOSFETS. 

5.1 SUGGESTIONS FDR FURTHER WORK: 

a) Initially the device chosen was 74HC04, The I-V 
characteristics (ref er Fig.is) of this high speed CMOS were 
found very different from that of normal CMOS 74C04, CD4007. 
Meyer model was found wanting in explaining the phenomenon 
observed^ A search for a model that can describe the I-V 
characteristics of 74HC04 can be carried out, and its 
parameter extracted. The I-V data for 74HC04 has been 
plotted in Fig. 13 for ready reference, 

b) The device under study does not provide for studying 
the effect of source to substrate bias. Since Meyer's modiel 
explains the effect of source to substrate bias, those N 
and P channel device (MFE 3001, MFE 3003 respectively) could 
be chosen which facilitate the application of source to 
substrate bias, and its behaviour studied. 

c) Observations were also made to determine the values 

of using a C-V plotter. It was found that the 
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transitio^v capacitance of the pro-tective zenner coinpletely 
masked the vaiues o-f CGS and CGD , Meyex's model of 
voltage dependent capacitances finds favour with latest simu- 
lation programme like SPICEl and SPICE2, hence the same may 
be analysed using N and P channel mosfet that have no zennex 
protection; as it would help throwup infoimation about 
fast surface state densities. The C-V data taken for N and 
P channel are shown in Figs, 14 and 15 respectively. The value 
for for N channel MOSFET expected is in the vicinity of 0.2 
pf and that of P channel MOSFET is in the vicinity of ,14 pf. 

d) In the noise measurement scheme data has been taken 

for only the strong inversion case. The weak inversion case 

involves measurement of subthreshold current measurements 

which are of the order of nA and less. An electrometer may 

2 

be employed to get for weak inversion region* 

There are many more facets of the device that ono 
could be interested in, but it is hoped that the above 
suggestions if implemented may just about roundup the 
d,c, and noise behaviour of the device under study. 
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APPENDIX A 

DEVICE PHYSICS CONCEPT REVISITED 

Th€ most important paraiBeter in semiconductor 
physics is the band gap E which separates the lowest energy 
state in the conduction band from the highest energy state in the 
valence band* 

For a intrinsic semiconductor the number of occupied 
conduction band level is given by 

n = ; N(E)F(E)dE (1) 

Ec 

where 

_ energy at the bottom of conduction band 

£ a: energy at the top of conduction band 
top 

N(E) =* density of states 

p^£^ _ Fermi-Dirac distribution function given by 


F(E) * 

l-fexp(*^^) 

where 

K = Boltzman’ $ constant 
T =s absolute temperature 


Ep = Fermi potential 



S Ml 


pm avMV'Cvt/kTM 
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From equations, (l) and (2) we have 


■>1 “ 

c 


txp (• 



(3) 


where is the effective density of states in the 
conduction band. 

Similarly the hole density near the top of the valence band 
is given by 


P 


Ny escp { 



(4) 


Also, the Fermi level Ep fox the intrinsic semiconductor is 


given by 


E -fE N , 

C- -V + is. In 
~ ^ 2 


(5) 


Hence the Fermilevel Ep for the intrinsic semiconductor lies 
very close to the middle of the band gap. Refer FigJ.6(a). 

The intrinsic carrier density is obtained by 

n.p = n^^ =» exp (-j^ ) (6) 

when impurity are introduced, the Fermi-level must adju 
itself to preserve charge neutrality. Consider the case 
Shown in Fig.l6(b) where donor impurities with concentration 
are added to the crystal. To preserve the electrical 
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neutrality the total negative charges (electron and ionised 
acceptors) must equal the total positive charges (holes and 
ionised donors) or 

n « •h p (7) 

where 

n =5 electron density in the conduction band 
p « holes density in the valence band 
Np * the no. of ionised donors 



Nq C 1 - 



( 8 ) 


Similarly when acceptor impurities of concentration 
N^( cm ) are added 



where 

g M ground state degeneracy of donor/acceptor 
g » 2 for donor impurity in ( and S^) 

* 4 for accepter impurity in (Y^ and S.) 

when impurity atoms are added, the n.p# product is still 
given by eqn. (6). At relatively elevated temperature most 
donor and acceptors are ionised. 
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n + » p 


(10) 

From equation (6) 

and (10) concentration of 

electrons 

holes in a n type 

semiconductor 


"no. - i C(VNa)+«Nd-Na)^ + ■* "i'] 

(11) 

« 2 

2 


*^i 

11 . 


^no. n^^ 
no 

Nd 


and Eq— E jj « KT In 

(-^) 

(12) 

Ep-E^ « KT In ) 

(13) 


Similarly the concentration of holes and electrons in a p 
type semiconductor is given by 


Ppo “ i f(N^-l^*4r,.2) 

(14) 

^2 „ 2 


P® Ppo 

(15) 

Ep*^Ey * KT In (j;j ) 

A 

(16) 

p 

Ej^-Ep « KT In (-^) 

i 

(17) 


Consider a metal insulator semiconductor structure of Fig. 
The thickness of the insulator is d and v is the applied 
voltage on the metal plate. 


IT. 
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Let 

j ^ » surface potential at the semiconductor surface 
w » r # t ♦ 

^ g *> potential difference between the Fermilevel 
Ep and intrinsic Fermilevel 

when th« MIS structure of Fig .17 is biased with positive or 
negative voltage, basically three cases may exist 


i) Accumulation 
ii) Depletion 
tli) Inversion 

The effect of applied bias and consequent band 
bending is illustrated in Fig. 4. Briefly we have the 
following cases fox a p type semiconductor 


IS 


h 8 


< 0 accumulation of holes 

= 0 flat band condition 


> >0 depletion of holes 

SB '^'i's 

(band bend downward) 
mid gap with n^ = Ps^^^i (i^transic 
c 0 nc entia t ion) 

At th. .utf.ce th. ewotron and hole density are given by 




'h'B 


Ps * Ppo 


(13) 

( 19 ) 
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B « ^ 

q 

n ^ 3 equilibrium electron density 

po 

• equilibrium hole density 
po 

is essentially the measure of bend in at 
the surface of the semiconductor 

Figa.7 shows the band diagram of a MIS structure with 
the band bending of the semiconductor (p type). The charge 
distribution is shovm in Fig,lg(b), For charge neutrality 
of the system it is required that 

- Qn + q NgW = (20) 


Where Qw * chaxge/unit area on the metal 

=s felectron/unit area in the inversion layer 
q N W » ionised acceptor/unit area in the space 
charge space charge width W and 
Q « total charge/unit area in the semiconductor 

In absence of any work function difference the applied 
voltage partly appear across the insulater and partly 
across the silicon. Thus 


and 


V w g 

I Sll- 

Vi •= t-id * ' £ i 'ti 


( 21 ) 

( 22 ) 
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Another quantity of interest is the turn on voltage at 
twhich the strong inversion occurs and is given by 

Vy( strong inversion) =» ^ + 2 'g (23) 

Also for strong inversion 

u' 8 T ^ y p C24) 

Also at the on set of strong inversion 
Qjj ** q N^W from eqn. (2o) 

The turn voltage therefore 



also since the surface potential 



the space charge width 


( 26 ) 



The above discussion helps provide the physics of MOS 
device and provides the basis for formulation of any 


model* 
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APPENDIX B 

THE LEVENBEEG-MARQUARDT ALGORITHM 

The* aigorithm ossontially is the implementation of matrix 
aquation 

(J^J Ad^)Aa « / [Y-Y*J 

STEP DESCRIPTION 

0 Enter bEiJi i»l,2>,.,n, the starting values of 
tho parameters. 

Enter NN, the no. of data points. 

Enter minm, a small no. such as the machine 
precision to which the Marquardt parameter# lambda#, 
is to be set# should it have underflowed. 

Let inc * 10 be the factor to increase lambda 
Let doc » 0.4 be the factor to decrease it 
Let lambda » .001 be the initial value of 
Marquardt parameter. 

Let phi * 1 be th® value of quantity used by NASH to 
avoid singularity. 

Enter EPS =» The convergence criteria for least square 
objective function. 

Enter GAMA * The convergence criteria fox change in 
parameters between successive iterations. 
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SXEi* DESCRIPTJOM 

Enttr function 

Ent«r analytical darivativaa 

1 Computa S tht currant objective function 

2 Ccmputa j^j, 

3 Coraputa [Y-Y*] thru 'function* 

4 Computa [Y-Y*]. 

5 Comput* {J^J + p^) 

6 Solve for A 

7 Compute ^ new * A old + ^ 

8 Compute R the new objective function. 

9 If convergence criteria satisfied ‘sstbp^' [l2]« 

10 If (R.LT.S) » dec and Go to 1, 

11 If (R.GT.S) «: inc and Go to 1. 

STOP 
END 


12 

13 



78 


TART 


SfHifr ttirtiMHrfd va(t»«« for N mkI P cMmoot 
Smfff TrMt MiiMtiorw for N wt^ P chonnot 

Soodfir SifurtNon oiwfkm for N mdP thomol 
.V«»f. MiMfloni for N sN P dmwol 


Appx C 


MfMlot 


VO-VW,* 

VooVDO 



VWVIt.tS 


V! VDO 



I — 

Mtmkm 

vaa-vhvse 

? CilMMlt 


Vi Vitf ti 


VOoVDO 


aofafxa MtaratioD i^n 
for Mdimnot, fM fiS* 


aR7arca 7i 
i>M SC 


/s| Vfl 


Miforc* triado Mpi for 
Mdianmi/M DS 


Mforu ui wv 
?MSO 







VoaVoogltS 
V ottx 








Sres 

V0>V»>^9 Mtta 









HE fMO?i CONFIGURATION (CO-A001) 


— 


j/p 


3,10.6 


h 


VOO 
H,2.11 
S 

P Channel 


I — 


0/P {1,5),12.(8.13) 


D 


h- 

I N Channel 

S 


rSG -19 b 


Pin CONFIGURATION CCO-4007i 





80 


APPENDIX D 

SPCCIFICATIONS ^ DEVICE UNDER STUDY 

The device under study is CD4007 (Dual complementary 
peir plus a CMOS inverter The pin diagram of the chip is 
as given in Fig« 20» 


The chip comprises of three complementary fsiirs 
of N channel and p channel enhancejoient mode MOSFETS* All 
inputs are protected by diode clamps For proper operation 
voltages at all pin must be constrained to (VSS-.3V) and 
(VDS ^ .SV). The supply voltage range can vary from 
3.0V-15V. 


The physical parameter of the device are listed 


below. 


N Channel 

W » 7*5 mil 

L =» 0*15 mil 


7^^ - 1200A^ 


420 c 


n?/V 
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tv E ^ o 

K = = .30«-3 A/'T 


sec 
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'OX 

^si 
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12 
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^ox 



For P _Cbannel 


W » 19.5 rail 


Vtx = 0-63 


0.25 mil 


Nr 


Chi = - 0.62/at T=300°K 


1200° A 
ox 


K = = .189E-3 \iA/yr 


2 

\x s= 165 cm yV sec 

N = 3X10^^ cm*"^ 


K » lax E N e 1.082 
P Eox 
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